On the theory of thermoelectric phenomena in superconductors 
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The theory of thermoelectric effects in superconductors 
is discussed in connection to the recent publication by Mari- 
nescu and Overhauser [Phys. Rev. B 55, 11637 (1997)]. We 
argue that the charge non-conservation arguments by Mari- 
nescu and Overhauser do not require any revision of the Boltz- 
mann transport equation in superconductors. We show that 
the charge current proportional to the gradient of the gap, | A 
found by Marinescu and Overhauser, is incompatible with the 
time-reversal symmetry, and conclude that their "electron- 
conserving transport theory" is invalid. Possible mechanisms 
responsible for the discrepancy between some experimental 
data and the theory by Galperin, Gurevich, and Kozub Pis'ma 
Zh. Eksp. Teor. Fiz. 17, 689 (73) [JETP Lett. 17, 476 
(1973)] are discussed. 



I. INTRODUCTION 

The purpose of the present paper is to discuss some as- 
pects of the kinetic approach to the thermoelectric pro 
erties of superconductors. As early as 1944 Ginzbur, 
suggested that in the presence of a temperature gradi- 
ent, there appears in a superconductor a normal current 
of the form given by 

j„ = -aVT. 

It was also pointed out by Ginzburg that the total cur- 
rent in the bulk of a homogeneous isotropic supercon- 
ductor vanishes because the normal current is offset by a 
supercurrent js so that the total current in the bulk 

jri + js =0. 

This makes impossible the direct observation of the 
thermoelectric effect in a simply-connected homogeneous 
isotropic superconductor. Ginzburg considered also 
simply-connected anisotropic or inhomogeneous super- 
conductors as systems where it is possible to observe ther- 
moelectric phenomena by measuring the magnetic field 
produced by a temperature gradient.,. 

As indicated in Refs. (see alsoQ), the best way to 
observe thermoelectric phenomena in superconductors, 
in particular, to measure the thermoelectric coefficient 
a is to make the superconductor a part of a bimetallic 
superconducting loop that may also contain weak links. 



Using the approach based on the Boltzmann equation 
for the normal excitations, the calculation of the coeffi- 
cient a for impurityj-scattering has been made in Refs. |^,^ 
- see also reviewg3ia. The expression for a has been later 
rederived in Ref. ^ using the Green's function version of 
the nonequilibrium statistical operator approach. In that 
paper the role of paramagnetic impurities was also dis- 
cussed. Based on the same method, an enhancement 
of the thermoelectric flux in superconductors contain- 
ing nonmagnetic impurities with localized states near the 
Fermi energy was predicted in Ref. |l^, see also Ref. |l|. 
As was proved in Refs. [9|-pd|, the expressions for a, ob- 
tained in Refs. |^,^ for the case of nonmagnetic impurities 
remain valid for an arbitrary relation between the coher- 
ence length of the superconductor and the electron mean 
free path. ■— ■ 

Recently, Marinescu and Overhauserll3 have proposed 
another method to calculate the transport coefficients 
a. In their approach, the principal contribution to the 
thermoelectric effect in superconductors comes from the 
dependence of the superconducting gap A on the tem- 
perature. For some typical interval of temperatures and 
impurity concentrations their results differ from that of 
Refs. by several orders of magnitude. Therefore, it 
is desirable to discuss the validity of their results. In 
the present paper we compare these approaches. We also 
briefly discuss how the theoretical results are related to 
the existing experimental data. 



II. THEORY OF MARINESCU AND 
OVERHAUSER 

Marinescu and Overhauser in Ref. ^ have proposed 
a method which they call "electron-conserving transport 
equation" . They introduce distribution functions, ^kt 
and 5-ki which differ from the distribution functions 
for the BCS excitations /p|,x (below, the spin index is 
dropped). 

The nonequilibrium part of the distribution function 
[see Eq. (47) of Ref. |l| is 
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k-VT. (1) 



Here, as in Ref. |T|, /3 = l/fcsT, Eu = y^A^{T) + el, 
and fk = (^ef^Ek + ^ ^hile efe = h'^P/2m -epis the 
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one-electron energy measured with respect to the Fermi 
level, ep- 

The relaxation time Ts is related to the relaxation time 
T„ for impurity scattering in the normal state: For quasi- 
particle transitions from k to k', 

= r^^|£^/c/efc|(wfcUfc' - WfcWfc')^ • (2) 
For future convenience, we write Eq. (p|) as 



= — ^ ( — I k • VT . 



(3) 
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Accordingly, the electric current density is split as j = 
j(-f) + j(-f^) . where 



with 



a 



HI) _ 4ejV(0) 



dcTsie + ep) 



keT^fA dA 
with iV(0) being the density of states per spin. 

III. DISCUSSION 

Now we are in position to discuss the results of Ref. 
The first term, Sg^^^ in Eq. comes from the coor- 
dinate dependence of the temperature T entering the 
Fermi-Dirac distribution function, and therefore, is of 
true nonequilibrium origin. A nonequilibrium term exists 
in the distribution, function found from the Boltzmann 
equation approachnfl (cited as Eq. (21) in Ref. |lj) with 



a very important difference: 



is even under e 



-e so 

that electrons above and below the Fermi surface do not 
tend to compensate each other as in Refs. |,|. The oppo- 
site symmetry in e is the reason why the thermoelectric 
current obtained by Marinescu and Overhauser is some 
five orders of magnitude larger than that in Refs. 0,0. 
The origin of the second term, ffjj. , is the r- 

dependence of A = A[(T(r)] as is obvious when 5g^^^ 
is identically written in the following form: 



Tsfk 



A(r) 

El 



d_ 
dr 



A{r] 



hk 

m 
(II) 



This means that in the approach of Ref. |12|, 5g\^ ' would 
exist irrespective of the origin of dependence of A(r) on 
the coordinate r. For instance, such dependence may be 



due to the variation of the chemical composition of the 
superconductor or to the spatial variation of the strain. 

Even if the dependence A(r) is due to one of these 
equilibrium mechanisms, the theoryll3 nevertheless pre- 
dicts the current 



or 



(4) 



where [3 = a^^^> (dA/dT)-^ . In our opinion, such a cur- 
rent is forbidden. Below we give physical considerations 
supporting this statement. 

It is well-known that any linear response process can 
be classified as either reversible or irreversible (dissipa- 
tive)E3 (see Ref. ^ for a general discussion of nonequi- 
librium thermodynamics of superconductors, including 
thermoelectric phenomena). The basis for the classifi- 
cation is the time reversal symmetry (T-synimetry). For 
instance, the charge current j induced by the electric field 
E, in a normal conductor is irreversible. Indeed, the 
current changes its sign under time reversal whereas E 
remains intact, and the Ohm law, j = crE, is not invari- 
ant relative to the T-symmetry transformation. Ther- 
moelectric current, j*^"^^ — —a'VT, is irreversible because 
the left-hand side changes its sign under time reversal 
whereas VT does not. Another example is the super- 
current, js — 2eNs\'s, proportional to the density of the 
Cooper pairs, Ng, and their velocity, Vg- This relation 
is T-invariant since both current and velocity are T-odd 
quantities. Consequently, supercurrent is reversible and 
compatible with thermal equilibrium. 

From this point of view, the current in Eq. (|^) is ir- 
reversible. Indeed, the gap function A {— |A|) is un- 
changed by the time reversal transformation ijj ip*, 
whereas the current changes its sign. Unlike the super- 
current, the irreversible current in Eq. (^ if existed would 
be accompanied by a steady entropy production. Being 
incompatible with equilibrium, the current must be equal 
to zero. 

The contradiction with the T-symmetry arguments can 
be also demonstrated by the following gedanken experi- 
ment. Consider a ring built of two superconducting arms, 
left and right, both in thermal equilibrium. The arms 
are thick enough, so in their bulks the magnetic field is 
completely screened. The right arm is made of a chem- 
ically inhomogeneous superconductor with a position- 
dependent gap, A^'''(r), varying along the arm from Ai 
to A2 if one moves counterclockwise. Fig. |^. The left arm 
is built of a homogeneous superconductor with the gap 
A^'^ The expression for the electric current, according 
to Eq. (0) and the theory of Ref. 12, reads 



J = 
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-/JVAW, right arm, 
left arm. 
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FIG. 1. On the equilibrium non-quantized magnetic fiux. 

Thus, there is an additional current in the right arm 
given by Eq. (||). Since the current vanishes in the bulk 
of the ring, i = 0, from Eq. (|) one gets the following 
equation for the phase X: 

(27r)-iVx^'''^ = ^o^^A + (/3m/2enA^M)VAW(r) , 
(2^)-iVx('^) =$o"^A. (6) 

Integrating the phase gradient along a contour C in the 
bulk of the ring shown in Fig. ^ and remembering that 
the phase gets an increment 27rn , one gets for the total 
flux in the ring, 

$ = n$o + ^^'^^o) , n = 0,±l,... (7) 

where the additional flux $(^^'^) is 

The lowest-energy state (at least for small values of 
(j)(*^0)^ corresponds to n = 0, and, therefore, the equilib- 
rium flux of the ring is predicted to be finite and equal 
to $(*^'^). A similar conclusion was made in Ref. 

Clearly, the finite magnetic flux and a finite electric 
current generating the flux arc not compatible with the 
underlying time-reversal symmetry of the Hamiltonian 
of the superconductors. It is true that the state of the 
system may have symmetry lower that the Hamiltonian: 
A ferromagnet gives an example of a system with spon- 
taneously broken T-symmetry. However, the underly- 
ing T-symmetry guarantees the existence of two macro- 
scopically allowed states related to each other by the 
T-transform: for a magnet, theses are states with re- 
versed magnetizations. In other words, the underlying 
T-reversal symmetry demands permission of any sign of 
the parameter that quantifies the spontaneous symmetry 
violation. 



In the case under conside|Cation, the sign of the flux 
in the ring is predeterminedtJ, given the geometry and 
the material parameters. This contradicts the time rever- 
sal symmetry, by which the existence of two equilibrium 
states with the opposite magnetic fluxes is compulsory. 
Therefore, the current in Eq. (Q), which is the source 
of the spurious flux, must be equal to zero. (Inciden- 
tally, a finite <I>(*^'^) would mean that there are material- 
dependent corrections to the flux quantization phenom- 
ena.) 

These are our general arguments for why we think that 
the current given by Eq. (m should not exist. Of course, 
the absence of the equilibrium current which would be 
solely due to a spatial dependence of |A| is well-known 
in the microscopic theory of superconductivit}il3. Tbeije^ 
fore, in our opinion, the predictions of the theor5il3^t3 
contradict the general principles, and this is why we be- 
lieve that the transport theory suggested in Ref. |l^ is 
invalid. 

Although the derivation of the "charge conserving 
transport equation" has not been presented by Mari- 
nescu and Overhauseili3 in enough detail, let us try to 
specify those part in their calculations which has lead to 
the above contradictions with general principles. In our 
opinion, it is mostly related to the transport Eq.(39) of 
their paper. 

The equation is formulated for the function (7kt which 
gives the occupation number of bare electrons [see 
Eqs. (37)-(38) of Ref. The transport Eq.(39) is iden- 
tical to that which would be obtained b% the standard 
procedure based on the Liouville theoremllj. As it is clear 
from Eq. (46), the Hamilton function H{r, p), which de- 
fines Tik — ~ (dH/dr) and f = {dH/dhk), is taken in 
Ref. |l^ to be the BCS excitation energy 

H{r,nk) = Ek{r) = ^'e^ -f |A(r)|2 . 

We believe that this procedure is qualitatively unsatisfac- 
tory: Due to the electron-hole quantum coherence, the 
bare electrons are not good semiclassical eigenstates, and 
the motion of wave packets built of them cannot be de- 
scribed by the Hamilton equation, even approximately. 

The main objection of Marinesao, and Overhauser to 
the Boltzmann equation approachaLl is based on an ap- 
parent non-conservation of "bare electrons" and thus of 
electric charge in the course of the propagation of a wave 
packet. As is well known, the BCS Hamiltonian itself 
does not support detailed conservation of the "bare elec- 
trons" while the charge conservation takes place only af- 
ter the quantum averaging over the quantum states pro- 
vided the complex pair potential A ~ |A|e'^ satisfies the 
self-consistency condition (see Appendix-for the details). 
It has been specifically emphasized inQ'El that the conti- 
nuity equation, 

dN , . 

e_+divj = 0, 

bears no explicit relation to the Boltzmann equation. In 
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superconductors, the charge conservation law plays the 
role of a subsidiary equation allowing one to find the 
phase X of the superconducting order parameter. Again, 
the conservation holds only as an average over different 
states and the self-consistency equation plays a crucial 
role to support the charge conservation. Thus, the objec- 
tions of Marinescu and Overhauser do not undermine the 
Boltzmann equation approach. Furthermore, the same 
result for the transport coefficient was derived using the 
Green's function method.El 



IV. EXPERIMENTAL SITUATION 

Let us briefly discuss experimental results. The ther- 
moelectric flux through a closed loop has been first mea- 
sured by ZavaritskiO. The results exhibited no serious 
discrepancies with th^hepry of Refs. ^,||. However, in 
the later experimentsE3 Ell the observed thermoelectric 
flux was much larger than the estimates from Refs. |^,^. 

In our opinion, the notable discrepancy between dif- 
ferent experimental results as wellpjpjejt*\pen some of those 
results and the microscopic theorjfl'El'Bil^l remains a chal- 
lenging problem. We present a brief review of different 
suggestions regarding this problem, see also p. 1856 of 
Ref. |2| where a discussion of possible reasons for the dis- 
crepancy is presented. 

Two possible suggestions can be made regarding the 
source of puzzling discrepaiicies between the theory and 
the experimental resultalacd. The first is focused on 
the differences between the realistic circuits, used in 
Refs. |l8|-^ and the simple theoretical modelo'El. 

The complications can arise, in particular, from the 
near-contact regions. In Ref. a thermoelectric loop 
consisting of an impure branch with higher Tc and a pure 
branch with lower Tc (passive and active branches, re- 
spectively) was considered. It was shown that if (i) the 
electronic thermal conductance of the passive branch is 
much smaller than that of the active one, and (ii) there 
is finite thermal flux through the contact between these 
branches, then there exists a large contact thermoelec- 
tric contribution to the measured flux due to the phonon 
drag. The reason is that the phonon thermal flux in one 
of the materials can not be transformed to the electronic 
one abruptly at the contact. Rather, the transformation 
takes place within a near-contact layer of finite thickness 
where the phonon flux within the active branch is much 
larger than that in its bulk. As a result, the contact 
contribution can exceed the predictions of Refs. ||J| by 
a factor ep/Qo where is the Fermi energy while 
Qjj is the Debye energy. However this enhancement, al- 
though substantial, seems to be still too small to explain 
the magnitude of the effect reported in Refs. pO| , pT| . 

Another suggestion is that the effects observed in 
Refs. p"8|-^ can be related to some temperature- 
dependent magnetic fluxes produced by external sources. 
A possible effect of such a sort was suggested in Ref. pi 



and later considered in detail in Ref. |2j. This effect is re- 
lated to a spatial redistribution of a background magnetic 
flux due to the temperature dependence of the London 
penetration depth. A, of the superconductor. Since the 
redistribution effect is proportional to \/L, where L is 
a size of the circuit, while the "true" thermoelectric flux 
within the superconducting circuit is proportional to a 
smaller factor, (A/L)^, even very weak background fields 
can produce a temperature-dependent flux. Indeed, the 
effective near-surface area ~ \L becomes rather large 
in the vicinity of where A > I /im. Then, the non- 
screened magnetic field of the Earth, for instance, may 
generate a temperature-dependent magnetic fiux as big 
as ~ 103$o- 

It is worthwhile to note that if the diameter d of the 
wires forming the loop is much less than L the "redistri- 
bution" effect can be suppressed by a small factor d/ L. 
The reason is that the contributions of the "inner" and 
the "outer" parts of the wire to the redistribution effect 
have opposite signs. We believe that it is because of the 
suppression of the redistribution eflfect that the first ob- 
servation of the thermoelectric fluxliH exhibited no serious 
discrepancies with the theory. 

To avoid the "redistribution" effect the authors of 
Refs. I^^exploited the experimental set where the ther- 
moelectric circuit had a shape of a hollow toroid, the ther- 
moelectric flux being concentrated within its cavity. The 
measuring coil was winded along the toroid. The main 
idea was that the background fields were screened out by 
the bulk of the toroid. In our opinion, there still existed a 
region between the coil and the bulk of the toroid where 
the background field can penetrate. In particular, this re- 
gion included the near-surface layer of the thickness ~ A. 
Thus the measured flux obviously included a contribu- 
tion of this " outer" region. Consequently, a temperature 
dependence of the penetration length brought about a 
temperature-dependent contribution provided that there 
existed SQpp background field. The fact that, as reported 
in Refs.EElO, the effect vanished for a homogeneous sam- 
ple, in our opinion, cannot serve as proof of the absence 
of the redistribution effect. Indeed, if the measuring cir- 
cuit is fully symmetric the contributions due to different 
parts of the toroid are compensated. However, for an in- 
homogeneous sample this symmetry will be absent, and 
the effect of the background field will be restored. 

We believe that to check,|experimentally the validity 
of the theoretical approachoQ it is practical, along with 
further studies of the thermoelectric ffux under different 
geometries, to study thermoelectric effects of other type. 
Among these effects, there is a specific interplay of a tem- 
perature gradient and a supercurrent in a superconduct- 
ing film. Due to such an interplay a difference between 
the populations of the electron- like and hole-iLkjaiffaaches 
of the quasiparticle spectrum is establishecOO'EltZl. As 
a result, a difference Ut between the electro-chemical 
potential and the partial chemical potential of the quasi- 
papticles appears. According to the experimental stud- 
ie^'c3, the measured values of Ut agree with the theo- 
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retical predictions. We would like to emphasize that the 
measurements of Ut are local and consequently are much 
less sensitive to the above-mentioned masking redistri- 
bution effect. Another way of local measurement is the 
thermoelectric modification of the Josephson effect in the 
SNS Josephson junction predicted in Ref. The theo- 
retical predictions obtained by the Boltzmann. equation 
approach agree with the experimental resultsEj. 
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V. CONCLUSIONS 



APPENDIX A: CHARGE CONSERVATION 



a 



In our opinion, the theory of thermoelectric effect in 
superconductors suggested by Marinescu and Overhauser 
in Ref. |l2| is not valid. In contradiction with general prin- 
ciples, the main contribution to the thermoelectric coef- 
ficient a calculated in Ref. |2| can be attributed to the 
spatial dependence of the order parameter A(r), which is 
due to the temperature dependence of the gap A(T) when 
T = T{r). As such, the spatial inhomogeneity of A does 
not disturb the thermal equilibrium. We have shown that 
the existence of an equilibrium current oc VA contradicts 
the time reversal symmetry, and, therefore, it cannot 
contribute to the thermoelectric coefficieBt. In addition, 
being irreversible, the current found inlia would lead to 
a steady entropy production in an inhomogeneous jS6 
lihrium superconductor, violating thermodynamicsEJ 
It is our belief that "the electron-conserving transport 
equation" proposed in Ref. |l^, from which the above un- 
physical results are inferred, is erroneous. nnn 

We have presented the conventional point of view,l3l3'Eil 
on the issue of the charge conservation in superconduc- 
tors: It is an intrinsic feature of the BCS mean-field 
theory that the charge conservation may be violated 
i.e. divj'n ^ 0, for any individual quasiparticle state 
Tpn- However, the total electric current j, which is the 
sum over the quasiparticle states is locally conserved, i.e. 
divj = 0, if and when the pair potential is self-consistent 
(see the Appendix for details). Therefore, the Boltzmann 
equation supports the local charge conservation, and we 
disagree with llae opposite claims in Ref. |l2|. 

The resultsoa of the Boltzmann equation approach 
are consistent with the experimental studies of the 
thermoelectrically-induced branch imbalances and cor- 
rections to the critical current of the Josephson junc- 
tionEj. The results of the thermoelectrically-induced 
magnetic fiux through a close topp exhibit substan- 
tial scatter,,-. Whereas the resultstll are consistent with 
the theorygq, the giant thermoelectric flux observed in 
Refs. |l8|-pi| is still not understood. We think that it is a 
challenging problem which is still open and it may require 
an account of additional sources of thermoelectrically- 
induced magnetic flux. However, we are convinced that 
within the framework of the physical picture involving 
quasiparticle diffusion, the Boltzmann equation approach 
of Ref. H does not in principle require any revision. 



A quasiparticle in a superconductor is described by a 
two-component wave function ijj — (") , with u and v be- 
ing the electron and hole components of the quasiparticle, 
respectively. The stationary wave function corresponding 
to the enesgy E is found from the Bogoliubov - de Gennes 
equationllj, 7i'0 = Exp, with the matrix Hamiltonian 



n = 



e(p - f A) 

A* 



C/(r) 



A 

^A) 



C/(r) 



(AI) 



Here ^(p) = p^/2to — ep For definiteness, we consider 
an isotropic s-wave superconductor so that the order pa- 
rameter A(r) does not depend on the momentum. Also, 
A is the vector potential, and U is the potential energy, 
e. g., due to impurities or the scalar electric potential. 

Generally, the eigenenergy E in Eq. (Al) may be pos- 
itive or negative. In the ground state (the condensate), 
the states with the negative energy are filled, and the 
positive-energy states are empty. As usual, the excitation 
is defined relative to the ground state, i. e. it occupies 
an ii^ > state or empties an i? < state. It is a prop- 



erty of Eq. (Al) that the eigenfunctions correspoftding 
the energy E and —E are related to each other agl3 





(-::) 


j and 







This property allows one to express the contribution 
of the negative energy states via the positive energy 
ones, and, thpiefore, exclude the negative energies from 
considerationll§E3. 
We denote 



Un(r) 

■y«(r) 



as the wave function of the excitation with the energy 
En > and fn,a as the distribution function of exci- 
tation {a =1,i being the spin); here n stands for the 
quantum numbers other than spiH-j^he observables can 
be expressed via u„, «„, and /„^CT-tJ't3 

The densities of charge, Q{r), and electric current, j(r), 
are given by the following expressions 



B„>0 

j= Yl (-l + /"T + /ni)j« (A2) 



-E„>0 
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where the partial charge, gn(r), and current, j„(r), den- 
sities are 



iigp) 



Re (u* vun 



is a conserved quantity i.e. divjn'^'''' = 0, for any solution 
to the Bogoliubov - de Gennes equation. The conserva- 

iqp) 



^-1; _ ^,*,->„ ^ -Tr— /^A \\ tion of j„ Meads e.g. to the conservation of probabilities 

m V c / m the Andreev reflection problem. 



Indeed, as it is rightly stated in Ref. ^ the effective 
charge of the excitation g„ is, generally, a function of the 
coordinate r. Accordingly, the total charge of a wave 
packet built as a superposition of V'n's varies while the 
packet propagates in an inhomogeneous superconductor. 
Of course, this violates the charge conservation on the 
level of an individual excitation. However, it is well- 
known that the local charge conservation is restored: (i) 
after summation over the states, i. e., one should con- 
sider only the total charge and current rather than q„ 
and jn, (ii) the pair potential A = |A|e*^ is taken self- 
consistently rather than as an input. Below, we elaborate 
upon this point. 

Unlike the exact Hamiltonian of interacting particles, 
the BCS effective Hamiltonian does not commute with 
the particle number operator because of the presence of 
the anomalous average term ACp^cLpi + h. c. For this 
reason, the charge conservation is n ot an automatic prop- 
erty of BCS-theory. From Eq. (|Ai|), 



divjn = -4Im A*u„v* , 



(A3) 



so that lack of detailed current conservation, i. e. div j„ ^ 
0, is obvious. 

For the total ele ctric current density j, one gets from 
Eqs. (0) and @, 



divj = 4ImA*(r)F(r), 



where 



F(r) 



E„>0 



(1 - /n,T - /«,i)«n(rX(l")- 



As discussed in Ref. ^ the pair potential A (A*) serves 
as a source (sink) of the charge. Again, the charge conser- 
vation is not guaranteed if the potentials in Eq. ( |AlD are 
arbitrary inputs. However, the Gor'kov self-consistency 
condition demands that 

A(r) = gF{r) 

where g is the coupling constant. If the complex potential 
A is self-consistent, one readily sees that divj = 0, i.e. 
the local charge conservation. 

These are our arguments against the point of view ex- 
pressed by Marinesku and Overhauser that the Boltz- 
mann equation scheme violates the local charge conser- 
vation. To avoid confusion, another point should be men- 
tioned. The lack of the detailed current conservation does 
not mean the absence of unitarity. Indeed, it generally 
follows from the Bogoliubov - de Gennes equation that 
the quasiparticle current j ^"^^^ , 
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